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ABSTRACT 


Cellular glaases are prime candi'^ste cMterials for the structural sub- 
strate of mirrored glass for solar concentrator reflecting panels. These 
materials possess properties desirable for this application such as high stiff- 
ness to veight ratio, dimensional stability, projected low cost in mass pro- 
duction and, importantly, a c^ose match in thermal expansion coefficient with 
that of the mirror glass. These siaterials are brittle, however, and suscepti- 
ble to mechanical failure from slow crack growth caused by a stress corrosion 
mechanism. 

This report details the results of one part of a program established to 
develop improved cellular glasses and to characcerize the behavior of these and 
commercially available materiils. Commercial and developmental cellular 
glasses were tested and analyzed using standard testing techniques and models 
developed from linear fracture mechanics. Two models describing the fracture 
behavior of these materials are developed. Slow crack growth behavior in 
cellular glass was found to be more complex than that encountered in dense 
glasses or ceramics. The crack velocity was found to be strongly dependent 
upon water vapor transport to the tip of the moving crack. The existence of a 
static fatigue limit was not conclusively established, however, it is specu- 
lated chat slow crack growth behavior in Region I may be slower, by orders of 
magnitude, than that found in dense glasses. 
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SECTION 1 


INTRODUCTION 


Cellular glasses have a long history as self-supporting insulation 
materials for both hot and cold applications. Therefore | optimizing their 
thermal insulating properties has been emphasized in developing these 
materials. Consequently, little interest was shown in the mechanical behavior 
of the cellular glasses because of the minima! structural requirements in this 
traditional application. 

Recent interest in achieving highly accurate reflective mirror systems 
for solar thermal energy conversion programs has resulted in the search for a 
durable, low-cost and light-weight structural substrate for mirror glass 
reflectors. Cellular glasses have been identified as a prime candidate for 
these optical, structural substrates due to their low density, close thermal 
expansion match to the mirror glass, and projected low cost. This novel 
application for these materials demands much greater structural integrity than 
the traditional insulation applications. An initial study was performed by 
JPL to assess the suitability of using this material in a mirror panel appli- 
cation (Reference 1). Subsequently, additional in-depth studies were carried 
out to gather basic design data and answer specific questions concerning the 
use of these materials in solar thermal systems. This report deals with the 
fracture mechanics investigations which include measurement of the critical 

stress intensity factor (K_.) for various cellular glasses and determination 

xc 

of the crack velocity vs. stress intensity factor relationship for Heavy Load 
Bearing (HLB) Foamglas® (trade name of the Pittsburgh Corning Corporation, 
Pittsburgh, PA.). Fast fracture strength determinations for this same material 
are also reported. 
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SECTION II 


THEORETICAL ANALYSIS 


Description of the stress intensity factor for dense, honogeneous, 
isotropic solids which exhibit linear elastic behavior can readily be 
described by the well-known relationship (Reference 2): 

Kf - or^ 

Where is the stress intensity factor for opening node I, <T^ is 

the stress applied to the material at some distance remote from the crack tip, 

Y is a dimensionless parameter which depends upon the geometry of the loading 

and crack configuration and C is the characteristic dimension of the flaw. In 

loading a material, once reaches a critical value denoted unstable 

crack growth occurs. K__ is typically considered to be an intrinsic material 

xc 

property which may be used to characterize the material's mechanical integrity. 

Cellular glasses are foamed, inorganic glasses (See Ref. 1 for a 
description of the manufacturing process). As such they are microscopically 
porous, heterogeneous materials on the scale of millimeters (Figure 1). In 
highly porous materials such as polymeric foams or inorganic glass foams, the 
physical interpretation of the terms in Equation I becomes uncertain. In such 
materials a microscopic structure exists which may profoundly influence the 
nature or interpretation of the variables in Equation 1. (The properties 
reported in this paper are bulk values.) The appropriate fracture mechanics 
model for cellular glasses which contain high volume fractions of pores 
(75-95 volZ porosity) may be one of modelling both the microscopic stresses in 
the cell walls and the nature of the cracks in these .:ell walls vis-a-vis 
modelling the cracks as though they existed in a continuous medium. This 
would be done with particular bulk elastic constants and nominal stress levels 
derived from the bulk elastic constants (see Figure 1). 

Figure la illustrates the bulk model for cellular glasses. This model 
considers the material as a composite material composed of two phases, one 
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Figure 1. Cellular Glass Models: Bulk and Microsco^ic 





being dense gleee end the other the gee in the closed poree. This aodel 
■uggeete use of the lew of aixturee to predict the Mchenicel properties. The 
pore gee occupies e specific frsction of the bulk voluae but does not possess 
sny Mchstixcsl strength or tensile elsstic aodulus. The bulk stress in this 
siodel is detersdned by dividing the spplied losd by the cross -sect ions 1 sres 
over which the losd is spplied, ss ■essured froa the bulk disMnsions of the 
speciaen. The bulk elsstic aodulus is then the bulk stress divided by the 
bulk strsin. The vsrious frscture aechsnics psrsaeters ere then deterained 
froa these bulk stresses end bulk elsstic constsnts. 

The aicroscopic aodel of celluler glssses is illustrsted in Figure lb. 

In this aodel the asterisl is considered ss s continuous, three diaensiousl 
structure coaposed of thin, doubly-curved, dense glees cell wells. In this 
aodel the stresses in the cell wells ere deterained by dividing the losd 
spplied to the body by the losd csrrying cross-sect ionsl sree of the cell 
wells. Tlte strsin in the cell wells is the rstio of the stress in the cell 
wells to the aodulus of the dense glsss in these wells. In this aicroscopic 
model.rotstions of the cell well structure, which produce bulk displsceaent 
under losd, asy occur without genersting spprecisble stresses in the cell 
vslls. If such flexing of the cell wells does occur, this %fould lesd to 
Isrger bulk displsceaents then would be predicted by the aicroscopic aodel snd 
one would expect devistions froa elsstic constsnts predicted by the Isw of 
aixtures. These devistions would csuse the aessured aodulus to be less then 
thst predicted by either aodel, especislly st lo%rer densities. 

Cellulsr glssses ere susceptible to the phenoaenon of subcriticsl creek 
extension (References 1 snd 3). As deterained for dense glsss, this is believed 
to be stress corrosion phenoaenon, the priasry corrosive aediua being wster 
vspor present in the environaent (References 4 snd S). In this physicsl pro- 
cess, when Che asterisl is subjected to s stress below its ultiaste strengtli, 
the flews inherent within the suterisl grow ss s function of tiae. The rste 
of growth of these subcriticsl flsws (i.e. the ersek velocity V) is s function 
of the stress intensity fsetor K^. Experiaencsl dsts on asny brittle cersaic 
asterisls hsve been sdequstely described by Che folUwing eapiricsl reistionship: 

V - AKj" (2) 


2-3 



This forn of the equation deacribing the V, relationahip tuffera froa the 
complication that the diaeneiona of A are dependent upon the atreaa eorroaion 
exponent n end the expreaaion ooea not contain the fundamental property of the 
material t A aiaple coapariaon of A valuea for auit<.^riala with different 

atreaa eorroaion exponenta ia not poaaible aince the phyaical baaia of the 
parametera A and n have not been eatabliahed and their behavior haa been 
obaerved to be coupled in experimental data. Likewiae, a compariaon of only n 
valuea betueen oMteriala ia not concluaive. Minnear and Bradt (Reference 6) 
and Sinea (Reference ?) have auggeated a more uaeful form of thia relationahip 
given by the following equation: 


V -/SiKj/Kjjj)** (3) 

In thia caae ia normalized to the value for the material and now 
haa the dimenaiona of velocity for all materiala, which greatly facilitatea 
compariaon of /9 and n valuea for different material a. 

The idealized relationahip of equation (3) ia illuatrated in Figure 2 

which ia plotted on a log va. log acale. Figure 2 indicatea a minimum required 

value of Kj for atreaa eorroaion to occur, i.e., a atreaa eorroaion limit 

below which no auberitieal eraek growth will oeeur. Above thia threahold 

value, K , the behavior haa been modeled (Refereneea 8 and 9) aa exhibiting 
o 

three regiona of eraek growth for denae glaaa and eeramiea. In region I, at 
low valuea of K^, the rate of atreaa aaaiated chemieal attaek at the crack 
tip controla crack propagation. Aa thia reaction ia atreaa aaaiated, increaa- 
ing the atreaa at the crack tip cauaea an increaae in crack growth velocity. 
Equation (1) illuatratea that thia proceaa ia autocatalytie under conatant or 
increaaing loading conditiona in that any crack growth (increaae in C) reaulta 
in an increaae in atreaa intenaity at the crack tip thereby cauaing an increaae 
in crack velocity. Aa the combination of increaaed crack length and/or 
increaaed atreaa cauaea an increaae in K^, the crack velocity increaaea in 
Region I, aa ahown by the curve, untii attaina the value at which the 
onaet oJ region II crack growth behavior occura. In region II the proceaa ia 
believed to be controlled by the rate of diffuaion of water vapor to the crack 
tip and crack velocity ia found to be nearly independent of atreaa intenaity. 
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Figure 2. The Variation of Crack Velocity as a Function cl 
Stress Intensity Factor for a Material Subjected 
to a Stress Corrosion Enviroment 


In region III, is close to 
sisting of theraally activated bond cleavage and again the crack velocity 
becomes very sharply dependent upon the stress intensity factor. 

For most engineering applications, the crack velocities of regions II 

and III are sufficiently rapid as to cause immediate failure of any component 

loaded to these levels. Region I behavior is, therefore, the region of 

interest for design purposes. There is controversy as to whether the threshold 

K does, indeed, exist or whether crack velocities continue to be finite at 
o 

lower K^'s but are merely of such a small magnitude as to be extremely diffi- 
cult to measure in the short time period available for experiments (Reference 
10), i.e., 10 m/sec. In dense glasses, has been determined by some 

investigators (References 11-14) to be on the order of 20Z to 40Z of the 
value for the material, that is, with very small decreases in K^, the crack 
velocity became negligibly small and was assumed to have a negligible effect 
on design lifetime. It is of engineering interest, therefore, to determine the 
static fatigue limit for the cellular glass material, since, at stress levels 
corresponding to Kj<K^ i.e. 0‘<cr^, designs are possible in which the 
static fatigue phenomena may be ignored in predicting reliability. This is a 
considerable simplification in the prediction of performance. 

This investigation attempted to determine both and for HLB 
Foamglas®, a soda-lime silicate cellular glass. Fast fracture behavior and 
slow crack growth behavior of the V vs. diagram depicted in Figure 2 were 
also investigated. Two other developmental cellular glasses were investigated 
to determine only. 


R.. and the process can be described aa con- 
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SECTION III 


EXPERIMENTAL PROCEDURE 


A. Kj^, DETERMINATION 

The critical stress intensity factor for opening mode, was 

determined for one coomiercial and two developmental cellular glass materials. 
The materials are listed in Table I which identifies the nominal chemistry, 
density and number of specimens available for testing. Single-edge notched 
beam (SENB) specimens, illustrated in Figure 3, were machined from large blocks 
of cellular glass received from the material suppliers. Each beam was weighed 
and dimensionally characterized for density determination and subsequent input 
to fracture mechanics equations. The beams were loaded to failure in uniform 
bending on a screw-type Instron Testing Machine (Model 1122, Instron Corpor- 
ation, Canton, Mass.). Four-point line loading was applied using 2.5-cm 
diameter cylinders. The specimen dimensions were chosen to preclude surface 
crushing of the cellular glass under the loading cylinders. 

The beams were loaded rapidly in an attempt to minimize the effect of 
slow crack growth during the loading interval. The Instron crosshead was set 
to run at 51 cm/min. Under this loading condition the interval of the fracture 
event from the initial load application to failure was on the order of 50 
milliseconds. Special instrumentation was required for accurate data 
acquisition due to the relatively slow response of the Instron chart recorder 
at these rapid loading rates. The signals from the load cell and from a 
linearly variable differential transducer used to measure crosshead motion, 
were input to a fast-response light-pen recorder (Honeywell Oscillographic 
Recorder, Visicorder Model 1858) set to run at 25 cm/sec. With this 
instrumentation the load data could be accurately acquired during the short 
time interval of the fracture event with a resolution of fractions of a 
millisecond. The critical stress intensity factor values were determined from 
the following equation (Reference 15): 
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Table 1. Nominal Oxide Composition^ Density and Number of 
Test Specimens of the Cellular Glass 


Experimental** 
Aliiminoboros il icate 
Cellular Glass 


Foams il® 75* 


Si 02 

72 

69 

67 

AI 2 O 3 

5 

11 

10 

B2O3 

— 

4 

4 

H 2 O 

— 

5 

5 

^®2®3 

— 

3 

3 

Alkali Oxides 
(R2O) 

13 

8 

11 

Alkaline Earth 
Oxides (RO) 

10 

— 

— 

Nominal Density 
Kg m“^ 

136 

320 

192 

Number of Specimens 

23 

13 

52 


*Trade name of Pittsburgh Corning Corporation, Pittsburgh, PA 

**Developmental material supplied by Pittsburgh Corning Corporation, 
Pittsburgh, PA. Coefficient of thermal expansion of 70 x 10"^/®C. 


IC 


^ V ,1/2 

bW^ 


.99 - 2.47(1) . 12.97(1)' 


23.17 




(4) 


where P^ is the critical applied load, is the four point loading moment arm 
length, b and W are specimen thickness and height respectively and C is the 
notch depth. the critical strain energy release rate can be calculated 

from the following relationship for plane strain conditions: 


IC (1 - 


. 2 ) 


IC 


( 5 > 
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2 

where E is the Young's modulus and u is Poisson's ratio. The (1- 1/ ) term 
is eliminated for conditions of plane stress. The difference between the two 
conditions is negligible due to the small magnitude of Poisson's ratio. 

Analysis of the dynamics effects during specimen loading and experimental 
verification of these conditions indicated that these dynamic effects con- 
tributed less than IX error in the worst case test. Acceleration forces on 
the top loading fixture and specimen were minimized by reducing the weight of 
the top fixture. The natural frequency of the fixturing and specimen was much 
greater than that required to strongly interact with the load application. 


B. DOUBLE TORSION TESTING 

A detailed analysis of the double torsion test technique may be found in 

the literature (References 9 and 16-22). Figure 4 illustrates the double 

torsion specimen as well as a definition of the various geometrical parameters 
and the four point loading configuration. The cell size of the cellular glass 
was on the order of 1-2 mm. Specimen sizes much larger than those reported in 

the literature as being useful for investigating dense glass and structural 

ceramics tiere utilized in this investigation in an attempt to achieve a repre- 
sentative microstructure in the test section. Two specimen conf igurations 
were used. The testing of thick section specimens was accomplished using the 
following typical dimensions: L ■ 46 cm« W * 30 cm^ * 13 cm, d ■ 5 cm 

and dw • 3.8 cm. It was found that to achieve an accurate compliance 
calibration a thinner section specimen configuration was required. In this 
configuration, L, W and remained the same while d « 1.9 cm and d^ * 1.3 cm. 

An elastic analysis of the torsional plate loading leads to the following 
relationship (see Reference 16): 


PW 


3 (1 ♦ i/) 
Wd^ d 


1/2 


( 6 ) 
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%rhere i/> Poisson's ratio and the other parameters are defined in Figure 4. 
Further analysis (Reference 20) has indicated that a correction factor *^(t) 
must be introduced due to the finite thickness of the double torsion specimens. 
The parameter t^ the thinness ratiO} is equal to 2d/W. The corrected equation 
then becomes: 


K 


I 



3 (1 ♦ i/) 

Wd^d f (t) 
w 


1/2 


(6a) 


Where ^(t) " 0.9212 for the thin section specimens used in this investigation. 
Equations 6 and 6a imply that is independent of crack length for this 
specimen configuration. Recent work (Reference 22) indicates that a constant 
independent of crack length, is obtained only when the crack is propa- 
gating in the central portion of the specimen i.e., end effects become impor- 
tant for relatively short or long cracks so care must be taken to achieve data 
free from these effects. 


The double torsion specimen may be tested in various modes to obtain the 
fracture mechanics data required for characterizing the subcritical crack 
growth behavior of a material. One mode utilized in this investigation was 
the load relaxation technique. If it can be demonstrated in a compliance 
calibration that the displacement, y, is linearly related to the crack length 
C and load P: 


y » P(BC ♦ X^) 


(7) 


where B and X are constants determined from the compliance calibration that 
o 

depend upon the elastic modulus of the material and the dimensions of the 
material and test device, then it can be shown that the crack velocity is 
related to the instantaneous values of the load and the load relaxation rate 
under fixed grip (constant displacement) mode of testing according to the 
following relation (see Reference 16): 


V 




( 8 ) 
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where the subscripts "1" and "f" refer respectively to the initial or final 
values of P and C for the test run. P and P are the instantaneous values of 
the load and the load relaxation rate respectively. The specimen is loaded 
rapidly to Pi where the machine crosshead is fixed. The crack begins to 
propagate rapidly but decelerates as the load relaxes according to equation 
(8). When the load relaxation rate decreases to sero the crack velocity is 
zero. In this test, one obtains data points at decreasing values of 
(K^/Kj^) on the V vs. (K^/K^^) curve of Figure 2 and approaches 
conditions i.e., V ■ 0, at long times. For the specimen configuration 
utilized in this investigation up to three load relaxation runs could be made 
on a single specimen to generate the plot in Figure 2. Once the load 
relaxation tests had been completed, the specimen was loaded rapidly to failure 
to determine P , i\e critical load from which K-_ could be determined from 

C Lw 

equation (6a). 
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SECTION IV 
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i 

I 
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RESULTS AND DISCUSSION 

A. CRITICAL STRESS INTENSITY FACTOR, 

The resulta of the single-edge notched beam tests are presented in 

Figure 5. The R-. data illustrate two standard deviations about the mean 
xc 

values of K-. And density as well as the number of beams tested for each of 
XC 

the three cellular glass sample populations tested. The solid curve is a 
least squares regression fit through the combined 96 individual ( Pf 
data points. The correlation coefficient for the linear fit is 0.981 which 
indicates an excellent fit to the data. 


The curve illustrates a definite linear dependence of measured bulk 
upon the specimen density over the limited range of densities investigated in 
this study. The material characterized in this program lies at the extreme low 
end of the density range; all densities being less than 15Z that of dense glass. 
It is of interest to note however, that extrapolation of the curve in Figure 5 
(see inset Figure 5) to the density of dense glass (2.55 gm cm ) results in a 


-3/2 


. This is well within the range of values presented 


value of 0.88 MNm 

by Wiederhorn (Reference 23) for of dense glasses measured in nitrogen gas 

o -3/2 

at 300 K. These values ranged from 0.749 0.011 tOim for soda-lime glass 

~3f2 

to 0.913+ .022 MNm for aluminosilicate glass. The macroscopic fracture 
toughness (bulk appears to be linearly dependent upon density for cellular 

glasses in the r«nge of densities studied. The dependence on density may 
be related to the for dense glass by the following law of mixtures rela- 
tionship between the two phases, dense glass and pore gas: 


(cellular glass) * (dense glass) 

DG 


(9) 
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(XIO MN m 



DENSITY (gm cm“^) 


Figure 5. Cellular Glass Fracture Toughness as a Function 
of Density 
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B 


DOUBLE TORSlONi V vt. (Rj/Rj^) 


Figure 6 ii the complience eelibration curve for the double torsion 
specinen. A point of special interest is that for dense glass and con- 
ventional ceramics «B so the X^/B tern in equation (8) becomes 
negligibly small. This is not the case for cellular glass of the specimen 
dimensions investigated in this study which were limited by the material 
supplied, X^/B was determined to be 16.0 cm which is the same order of 
magnitude as C. , ^ and therefore could not be neglected in the analysis. 

The slope of the compliance curve B, is related to the shear modulus / a by the 
following relation (see Reference 20): 


B 


3W 

m 

/xWd^^(t) 


( 10 ) 


From this equation /iwas calculated to be 4.43 x 10° Nm~^ (6.43 x 10^ psi). 
Under the assumption of linear elastic behavior, which appears valid in this 
case, Poisson's ratio i/, was calculated to be 0.17 from the following equation: 


1/ - L.1 > 


( 11 ) 


Where E is the sonically determined elastic modulus reported by the oianu- 
facturer to be i.o x 10^ x 10^ psi). 

In the analysis of the double torsion data the normalized values of Kj 
are of interest. From equation (6a): 



( 12 ) 


where Kj stress intensity factor at any load P, P^. is the critical 

load attained upon rapid loading to failure after completion of the load 
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COMPLIANCE (XIO”® mN’ 



COMPLIANCE 000“^ in. lb 





relaxation testa. From equations (1) and (12) one finds* 



TsTT' 


( 13 ) 


where o-^ is the distribution of the applied stress below which no slow crack 

growth w? '.1 occur and S is the distribution of the inert fast fracture 

strength, both being functions of the probability of failure of the oMterial, 

i.e., the flaw severity distribution in the naterial. Assuaing that the flaw 

distribution which detensines (S) is identical to that for (cr^), one suiy 

determine the static fatigue limit for the cellular glass from the 

and P values measured from the double torsion tests and from a weakest lirk 
c 

statistical analysis of the fast fracture behavior of the material. This 
distribution of (o^) would be the lower bound of region I slow crack growth 
behavior as a function of probability of failure. 

Seven thiwk section specimens (d ■ 5 cm, dw ■ 3.8 cm) were loaded to 
approximately 352>40X of where the crack propagated and was 
approached at the low end of the V vs. (K^/K^q) curve. *^e specimens were 
then rapidly loaded to determine P^. The ratio of (Kq/Kk;) determined 
from equation (13) ranged from 252 to 372 with an average of 312. This is in 
good agreement with the behavior exhibited by dense glass (see References 
11-14). 


A sample population of 49 unnotched speciswns of the coanercially 
available Foarnglsa^ soda-lime cellular glass was rapidly fractured in four 
point loading to determine the fast fracture failure stress d* stribution S. 

The specimen dimensions and loading configuration were as depicted in Figure 3 
with the exception that the saw^d notch was not present. Figure 7 is a plot 
of the two parasMter Vfeibull statistical analysis of the experimental data. 
Figure 8 is a plot of the probability of failure as a function of where 
from equation (13)t 


(or ). - 0.31 (S). 
0 1 1 


(14) 
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Figure 7. Foraglee® Feet Fracture Strength Weibull AnalyaU 
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PROBABILITY OF FAILURE 
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Figure 8. Foanglas® Weibull Probability of Failure vs. 
the Static Fatigue Stress Limit 
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Equations (8) and (12) were eaployed to analyse the alow crack growth 
behavior of the Foanglae® soda-liate cellular glass. Data analysis of the thick 
section specinens yielded invalid results due to the thickness of the speciisen. 
The double torsion analysis assumes a thin plate approaching a membrane. The 
thick specimens did not conform to this assumption. Thin section (d ■ 1.9 cm, 
dw ■ 1.3 cm) specimens with a correction factor for finite thickness of the 
specimen did result in valid data. Figures 9 and 10 present the slow crack 
gro«rth behavior of tvo such thin section specimens. Also included in the fig- 
ures are curves which are generated from the data of Wiederhorn et al. (Ref- 
erences 18 and 24) for dense soda-lime glass in various humidities, presented 

on the normalized axis of (K_/K _). In this normalization process a value 
5 -3/2 ^ 

of 7.5 X 10 Nm was assumed for of dense soda-lime glass. 

The three regions of slow crack growth depicted in Figure 2 are not 

apparent in the experimental curves of Figures 9 and 10. This may be a 

consequence of measuring only over a limited range of (R /K .). The 

I xc 

stress corrosion exponents (n). that is, the slopes of the curves shown in 
Figures 9 and 10, are plotted in Figure 11 in relation to the percent con- 
fidence interval computed from the experimental data. In a similar manner the 
intercept parameters, (/9), are plotted in Figure 12. The high values of n 
measured from these plots indicates a much stronger dependence of the crack 
velocity on (K /K ), i.e. steeper slope, than has been experimentally 
observed for dense glass except in Region III. Figure 12 indicates that the 
slopes of Figures 9 and 10 are indistinguishable above the 87t confidence 
intervals. This would indicate that there is some degree of variability in 
Foamgl as!® materials which is certainly plausible due to the inhomogeneity of 
the cellular glass structure. 

Diffusion controlled Region II behavior is not observed in Figures 9 and 
10. Perhaps this is due to the limited range of (R^/R^) over which the 
specimens were tested. However, several investigators (References 25-27) 
have noted that Region II behavior is a function of the specimen configuration 
and dimensions. Quackenbush et al (see Reference 27) have demonstrated that 
Region II behavior is comprised of a combination of environmentally dependent 
fracture in accessible regions near the sample edges which are accessible by 
diffusive transport of water vapor, and enviromsentally independent fracture 
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FiRure Q. Crack Velocity vs. Stress Intensity Factor for 
Foamglas'^ Spec imen 1 
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Figuro 10. Crack Velocity vs. Stress Intensity Factor for 
Foamglas*^^ Spec imen 2 
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Figure 11. The Stess Corrosion Exponent Intervals Determined From the Curves 
of Figures 9-10 for Various Confidence Levels 
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log STRESS CORROSION COEFFICIENT (m Sec 



occurring in the nore renote center of the fracture front. ThuSf as the 
specimen thickness decreases t the amount of material in the interior in which 
the crack is propagating via environmentally independent crack growth also 
decreases. In the limit of small thickness dimensions, the entire crack front 
becomes accessible to the corrosive environment and Region I behavior exists 
up to the point where the combined effects o' crack velocity and stress 
intensification are such that Region III behavior prevails. This process would 
eliminate Region II behavior altogether. This phenomenon could explain the 
absence of Region II in cellular glasses. The crack propagates as a thrrujth 
crack in the extremely thin (<0.1 mm) cell walls. Water vapor in the sur' 
rounding pore environment would always be accessible to the crack front, i.e. 
all diffusion paths from the aqueous cell gas environment to the crack front 
are short. The diffusion time required is small compared to the residence 
time of the crack front in the vicinity of the cell, thus environmentally 
assisted fracture behavior would exist up to velocities sufficient to intersect 
the Region III curve. Due to the thin cell walls, any appreciable Region II 
behavior %fould never occur. The slope of the defined Region I in Figure 10 is 
approximately 30. This is greater than the value of 16 usually measured for 
dense soda-lime glass but in close agreement with the value 31.2 
measured by Conley et al. (see Reference 3) on the sasie foamed glass siaterial 
using the dynamic fatigue technique in flexure tests. 

The magnitude of at the transition between Region III and Region I 
crack growth behavior may not be analogous to those seen in dense glasses. 

This transition may occur at considerably lower velocities than those observed 
in dense glass because of the reduced flux of water vapor transport to the 
cells surrounding the advancinr, crack. In closed cell cellular glasses, water 
transport from the surface ambient to the crack tip stust be through the cells. 
This is shown schematically in Figure 13. 

Effusion of water vapor from the ambient enviromaent through the small 
cracks in each cell wall is required. In addition, the water vapor may be 
readily combined with the resident cell gases (e.g., SO^) and adsorbed on the 
activated glass of the cell walls, effectively binding the water chemically 
and reducing the activity gradient for diffusion. The partial pressure of 
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Figure 13. Schematic of Water Transport to the Crack Tip Cellular Glass 


active water vapor in the cell atmoaphere would have to be built up in each 
cell to provide the driving force for diffusion of vapor into the next cell. 
This reduced rate of weter vapor transport oay result in a very low partial 
pressure of water vapor in the cells at the crack tip, even for cracks moving 
quite slowly. This dry condition at the tips of slowly moving cracks would 
result in a suppresion of Region I behavior at higher crack velocities %rhere 
it is nonsally observed in dense glass (see Figure 14). 

The knee in Figure 2 at K is of interest because it defines the lower 

o 

boundary for slow crack growth in siaterials. The data of Weiderhorn et al. 

(see Reference 24) indicate that the knee for soda-lime glass in water at 2S^C 

—8 —1 

occurs at crack velocities of about 10 m sec . Evans and Weiderhorn 

(see Reference 10) state that "...K occurs at such low velocities in ceramic 
-10 -1 ° 

materials, <10 m sec , that its existence has not been generally 

proven." If linear behavior is predicted from the range of (R^/Rj^) “ .31, 

determined from the low loading range previously discussed, extrapolation would 

yield crack velocities at the R knee in the curve on the order of from 
-15 —1 -40 -I ® 

10 m sec to 10 m sec . This has important ramifications for 

design purposes. If the design lifetime could allow Initial crack velocities 
-13 -1 

of 10 m sec for flaw growth of 3 mm, the curves would predict design 
(Ri/Ric) values of from 0.36 to 0.58 which would greatly increase the 
allowable design stresses for the material. 

The solution for R^^ from equation (6a) using the measured and U 

calculated from the compliance calibration yields a value of 5.85 x 10^ 

—3/2 

Nm ' . R^^ calculated from equation (6a) should be independent of the 

crack length of the specimen. R^^ calculated from equation (4) using single- 
edge notched bems data is strongly dependent upon crack length. The crack 
length used in equation (4) was the measured depth of the sawed notch. The 
machining process, by its very nature, introduces microflaws into the material 
thereby extending the length of the notch beyond the tseasured depth. The 
lengths of these extensions are unknown (possibly related to the cell size of 
the cellular glass) but one would expect R^^ values measured in notched beams 
to be some%dtat less than values measured in the crack length independent 
double torsion test. 
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Figure 14. A Comparison of the Slow Crack Growth Behavior 
for Cellular Glasses in Relation to that of 
Dense Glasses 



In any given number of 8BMB apecimena it is reasonable to expect that the 
l<;:igths of the microflaws produced ahead of the sawed notch by the machining 
process would be distributed over some range. In a very large number of 
specimens this sise range of small flaws ahead of the sawed notch would extend 
from sero to sosw maximum value. If is assumed to be a material 

constant) Equation 4 indicates that the effect of increaaing the crack length 
would be a decrease in the critical applied load. Since the experimental data 
analysis utilizes the measured critical load and the measured value for the 
sawed notch depth to determine from equation 4, the effect of an unknown 
increase in crack length due to machining microflaws would lead to a low 
estimate of the for the material. It is not important to determine this 
flaw size distribution but to realize that the notch depth measured on the 
specimen with the smallest microflaw ahead of it would be the closest approxi- 
mation to the actual notch depth of any specimen in that sample population. 
Testing of this specimen would result in the largest calculated value of 
and would therefore be the best approximation of the of the sample 
population. 

The maximum of measured on 29 Foamglaa^ SEME specimens was 5.57 x 

10^ Nm which is within 51 the value 5.85 x 10^ Dm measured from the 

double torsion specimen. The isaximum values measured on the two developmental 

4 -3/2 5 -3/2 

cellular glasses are 7.00 x 10 Nm and 1.20 x 10 Nm . The maximum 

values of measured from the three sample populations as a function of 

density are included on the plot of Figure 5. Extrapolation of the least 

squares fitted linear curve through these three points to the density value for 

-3/2 

dense glass yields a value of ■0.77 MNm which is in close agreement 
with accepted values of this quantity (see insert Figure 5). These maximum 
values of for cellular glasses, tend to lie at the high end of the 

density distribution for each sample population. Insufficient data existed to 
adequately consider the confounding effect of density variation within each 
population in this simplified extrapolation. An analysis which considers this 
density variation is expected to yield results similar to those presented here. 
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C. CRITICAL STRAIN ENERGY RELEASE RATE, 

Gj^ values were celeuleted from equation (5) ueing the maximum valuea 
of measured from each notched beam sample population as the best estimate 
for Young's modulus values were determined from Independent measure- 

ments of four point loaded unnotched beams. The calculated values of G_^ 

-2 

ranged from <!.2 to 4.8 Jm . Wlederhorn (Reference 23) reports values 
ranging from 7.6 to 9.5 Jm for dense glasses. The difference In these 
values Is unexplainable at present. 


D. CRITICAL FLAW SIZE DISTRIBUTION, C^ 

c 

The distribution of the Inherent critical flaw alee distribution for a 
sample population of cellular glass can be derived from equation (1) and a 
Welbull weakest link statistical analysis for the probability of failure as a 
function of fast fracture applied stress. This distribution Is given by the 
following relation: 



where C^ is the critical flaw size distribution, m is the Welbull slope, 6 la 
the scale paraaMter in the Welbull analysis, P is the probability of failure 
and the other paraawters are defined in equation (1). Three cases for Foamglas® 
specimens are analyzed, two on a euicroscale and one on a microscale. 

Case 1, on the aucroscale, assumes the bulk value estimate for R.^ ■ 

5.6 X 10^ Nm"^^^ (50.7 psl in m ■ 19.79 and ®«8.3i x 10^ Nn”* 

(120.42 psi) determined from the two parameter Welbull analysis of Foamglas® 
fast fracture tests, and Y ■ 1.128 for a penny-shaped crack (Reference 28). 

The resulting equation becomes: 


0.353 


*•» (rb) 


-O.IO 


cm 


(16a) 
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This c«t« would b« applicoblo to a half-ponny-thapod flaw on tha aurfaea of 
tha callular glaaa apaclnan. This distribution of critical flaw site ia 
illuatratad in curva 1 of Figura 15* 

Caao 2 ia also on a nacroacala and aaaumaa tha aana bulk valuaa for tha 
various paranatars with tha axcaption of tha gaonatric parasMtar Y. In thia 
casa Y • 1.965 corraaponding to an adga notch into a aemi-infinita plata 
aubjactad to tanaion (saa Rafaranca 2). The aquation for tha flaw sisa 
distribution in this case baconwsx 


C 


c 


0.114 



(16b) 


Thia ease approxinatea a flaw on the tensile surface traversing the width of 
the apacimen. It ia conceivable that a flaw conpoaad of connected open 
surface pores nay traverse the surface of the apacinan. Thia flaw size 
distribution ia illustrated in curve 2 of Figure 15. 


Case 3 treats the nodel, on a nicroscale, as a through surface crack in 

5 “3/2 

the thin wall of an individual cell. It asaunea K.. * 7.5 x 10 Nn 
1/2 

(682.4 psi in. ) for dense glass, n ■ 19.79 determined from the Weibull 

analysis and a value for 6 adjusted to account for the stress in the cell wall. 

A value for the stress in the cell wall is derived in Appendix A as being 27.8 

7 -2 

times the bulk stress value. Therefore 6 is estimated to be 2.31 x 10 Nn 
(3350 psi). The equation for this flaw size distribution becomes: 


C ■ 0.028 
c 



(16c) 


This distribution is illustrated in curve 3 of Figure 15. 

The three cases present models that are plausible. No attempt to verify 
or disprove any model has been attempted in this study.. The small magnitude 
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Glass Ual 








of the flaw size predicted by case 3, however, is subject to question. The 
machining process of sample preparation is expected to introduce flaws of this 
magnitude or greater. If the stress in the cell walls is actually higher than 
predicted by the model described in Appendix A due to "structural" effects in 
the microstructure, the flaw size distribution would be shifted to smaller 
sizes. 


The case 2 prediction is of a believable magnitude, on the order of the 
pore size. Conceivably the open pores could be connected across the thickness 
of the speciisen and v:uld lead to a through crack across the surface. Case 3 
predicts rather large ^enny-shaped flaws which have not been verified. 


SECTION V 


SUMMARY 


The bulk critical stress intensity factors of one coamercial and two de- 
velopmental cellular glasses were measured using single-edge notched beams in 
four points bending. A linear dependence of fracture toughness on density has 
been determined over the range of low densities investigated. Extrapolation 
of the linear curve to densities of dense glass shows good agreement with lit- 
erature values for of dense glasses. 

K f V vs. R. and R-_ determinations were made on the commercially 

O X xc 

available cellular glass utilizing the double torsion technique. R^/R^^ 
values ranging from 0.25 to 0.37 were observed to be in good agreement with 
observations on dense glass. The V vs. R^ relationship exhibited a much 
steeper dependence of the slope of V upon R^ than has usually been observed 
for dense glass. Region II diffusion controlled slow crack growth behavior was 
conspicuously absent indicating water vapor from the pore environment may al- 
ways be available at the tip of the growing crack. A physical model has been 
presented explaining this. The compliance calibration of the double torsion 
specimen provided a determination of Poisson's ration v- 0.17 and R^^ was 
determined to be 5.85 x 10^ Nm 

Microflaws are believed to be introduced by the machining process at the 
base of the notch in the SENB specimens used in this study. These microflaws 
result in the actual crack length being larger than the measured one by an 
unknown amount. This results in an incorrect, lower calculated value for R^^* 
The maximum calculated value for R^^ is believed to be obtained with the 
specimen which contains the mimimum size microflaw ahead of the machined notch 
and is therefore the best estimate for the of the sample population. 

The critical flaw size distribution has been analyzed for three models of 
the flaws in Foamglae^. Microstructural analyses have not been undertaken to 
verify which of the three cases is the active system. The model consisting of 
open surface pores linked across the width of the specimen appears to be more 
plausible explanation. 
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SECTION VI 


CONCLUSIONS 


This report details the findings of one part of a program to develop and 
characterize the behavior of cellular glasses to determine their suitability 
for utilization as a structural material in solar optical structures. Based 
upon the results of this study alone, a recommendation as to the durability or 
reliability of cellular glasses in optical structures cannot be made. Such a 
reconmendation must be made based upon the results of the entire program. The 
intent of the investigation reported herein was to make certain fundamental 
laboratory measurements on the materials. Because of time and funding limita- 
tions, these measurements were accomplished on only a small sample of mate- 
rials. No attempt was undertaken to evaluate the statistical variation of 
materials from different production runs. A follow-up study would be required 
to evaluate properties over months of production, to provide representative 
engineering design data on production materials. It is felt that the basic 
behavior and mechanisms determined in the present study should describe 
generally the behavior of any lot or batch of this material. 

The following may be concluded from the results of this study: 

o Slow crack growth behavior appears to be a more complex phenomenon 
in cellular glasses than in dense glass. 

o The existence of the stress corrosion threshold, K^, was not 
conclusively established for cellular glass. 

o Crack velocity at any Kj level may be dependent upon water vapor 
transport to the crack tip, more so than in dense glass or 
ceramics. 

o V vs. Kj crack growth behavior is dependent upon the crack 

lengths in the material due to the fact that water vapor transport 
is dependent upon crack length. It has not been necessary to 
consider this in dense glass and ceramics. 
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The data measured in this study provide a description of slow crack 
growth behavior in one regime of time, crack length, stress and resulting 
velocities. However, the behavior observed in this short study msy not be 
representative of the behavior of a cellular glass component in the field 
during a thirty year service life. To predict real service lifetimes with 
hi^ confidence, additional work will be required. Two options for increasing 
the confidence of prediction are available. The first is long term static 
fatigue testing. The second, requires a matrix of accelerated tests to 
determine the time dependence of water vapor diffusion to the crack tip. This 
testing would provide an understanding of the complex phenomena involved in 
the slow crack growth process which is necessary for predicting long term 
behavior. 

The identification of a stress corrosion threshold, K , was made 

o 

difficult by the fact that the material did not exhibit the typical stress 

corrosion behavior described as regions I, II and III. From the results of 

the present work, it appears that water assisted slow crack growth will be 

slower in the cellular glasses than in dense glass. If extrapolations to low 

crack growth velocities from the present data can be justified, design 

stresses on the order of 35% of the fast fracture stress (S), corresponding to 

“13 “1 

initial crack velocities on the order of 10 m sec , may be possible. 

It should be recognized, however, that the parameter (c/S) • 

not sufficient to determine the time required for a crack to reach a critical 

size. The time required for a crack to grow from initial to critical size is 

also dependent upon the initial crack length. Both the V vs. (K^/K^^) 

relationship and the initial flaw size distributions are required to predict 

lifetimes. 


6-2 


REFERENCES 


1. Giovan, M. , and M. Adams » Evaluation of Cellular Glasaes for Solar 
Mirror Panel Applications , JPL Publication 79-61. June IS, 1979. 

2. Paris, P.C., and G.C. Sih, ASTM Special Technical Publ. No. 381, 1965. 

3. Conley, P.H., Chandan, H.C. , and Bradt, R.C., "Dynamic Fatigue of 
Foamed Glass," Fracture Mechanics of Ceramics , Vol. IV , Edited by Bradt 
et al.. Plenum Press, New York, 1978. 

4. Charles, R.J., "Static Fatigue of Glass: I, II," J. Appl. Phys. Vol. 29 

No. 11, pp. 1549-1560, 1958. 

5. Wiederhorn, S.M., "A Chemical Interpretation of Static Fatigue," J. Am. 
Ceram. Soc., Vol. 55, No. 2, pp. 81-85, 1972. 

6. Hinnear, W.P., and Bradt, R.C., "(K-V) Diagrams for Ceramic Materials," 
ibid., Vol. 58, No. 7-8, pp. 345-46, 1975. 

7. Sines, G., "Rational! . Crack Growth and Time-to-Fracture of Brittle 
Materials," ibid., Vol. 59, No. 7*^, pp. 370-71, 1976. 

8. Wiederhorn, S.M., "Influence of Water Vapor on Crack Propagation in 
Soda-Lime Glass," ibid., Vol. SO, No. 8, pp. 407-414, 1967. 

9. Evans, A.G., "A Method for Evaluating the Time-Dependent Failure 
Characteristics of Brittle Materials - and its Application to 
Polycrystalline Alumina," J. Mat. Sci., No. 7, pp. 1137-1146, 1972. 

10. Evans, A.G. a.iJ Wiederhorn, S.M., "Proof Testing of Ceramic Materials - 
An Analytical Basis for Failure Prediction," Int. J. Fracture, Vol. 10, 
No. 3, pp. 379-392, 1974. 


7-1 



11. Holland^ A.J.« and Turner^ W.E.8.» "The Effect of Sustained Loading on 
the Breaking Strength of Sheet Glassy" J. Soc. Glass Teeh., Vol. 24, Mo. 
101, pp. 46-57, 1940. 

12. Baker, T.C., and Preston, F.W., "Fatigue of Glass under Static Loads," 

J. Appl. Phys., Vol. 17, pp. 179-178, 1946. 

13. Shand, E.B., "Experimental Study of Fracture of Glass: I, The Fracture 

Process," J. Amer. Ceram. Soc., Vol. 37, No. 2, pp. 52-60, 1954. 

14. Mould, R.E., "Strength and Static Fatigue of Abraded Glass Under 

Controlled Ambient Conditions: IV, Effect of Surrounding Medium," J. 

Amer. Ceram* Soc., Vol. 44, No. 10, pp. 481-491, 1966* 

15. Brown Jr., W.F., and Stravley, J.E., Plane Strain Crack Toughness Testing 
of High Strength Metallic Materials » ASTM Special Technical Publication 
Mo. 410, ASTM Philadelphia, Pennsylvania, 1966. 

16. Williams, D.P., and Evans, A.G., "A Simple Method for Studying Slow 
Crack Growth," Journal of Testing and Evaluation, Vol. I, No. 4, pp. 
264-270, July 1973. 

17. Evans, A.G., "A Simple Method for Evaluating Slow Crack Growth in 
Brittle Materials," Int. J. Fracture, Vol. 9, No. 3, pg* 267, 1973* 

18. Weiderhorn, S.M., "Subcritical Crack Growth in Ceramics," Fracture 
Mechanics of Ceramics, Vol. II , edited by Bradt et al.. Plenum Press, 

Mew York, 1974. 

19. Beachem, C.D., Kies, J.A., and Brown, B.F., "A Constant K Specimen for 
Stress Corrosion Cracking Tests," Materials Research and Standards, 

MTRSA, Vol. II, No. 4, pg. 30, 1969. 

20. Fuller E.R., Jr., "An Evaluation of Double Torsion Testing-Analysis," 
ASTM Special Technical Publication 678, "Fracture Mechanics Applied to 
Brittle Materials," pp. 3-18, 1979. 


7-2 



21. Plecka, B.J.t Fuller, E.R., Jr., and Koepke, B.G., "An Evaluation of 
Double Torsion Testing-Experimental," ibid., pp. 19-37, 1979. 

22. Shetty, D.R., and Virkar, A.V., "Determination of the Useful Range of 
Crack Lengths in Double Torsion Specimens," J. Am. Ceram. Soc., Vol. 61, 
No. 1-2, pg. 93, 1978. 

23. Wiederhorn, S.M., "Fracture Surface Energy of Glass," J. Am. Ceram. 

Soc., Vol. 52, No. 2, pp. 99-105, 1969. 

24. Wiederhorn, S.M., and Bolz, L.H., "Stress Corrosion and Static Fatigue 
of Glass," J. Amer. Cer. Soc., Vol. 53, No. 10, pp. 543-548, 1970. 

25. McKinney, K.R., "Effect of Specimen Size on Crack Propagation in Glass," 
Rep. NRL. Prog. No. 9, pp. 27-29, 1974. 

26. Chandan, H.C., Bradt, R.C., and Rindone, G.E., "Dynamic Fatigue of Float 
Glass," J. Am. Ceram. Soc., Vol. 61, No. 5-6, pp. 207-210, 1978. 

27. Quackenbush, C.L., and Frechette, V.D., "Crack-Front Curvature and Glass 
Slow Fracture," ibid., Vol. 61, No. 9-10, pp. 402-406, 1978. 

28. Sneddon, I.N., "The Distribution of Stress in the Neighborhood of a 
Crack in an Elastic Solid," Proc. Roy. Soc., London A187 , pg. 229, 1946. 


7-3 


APPENDIX A 


A MODEL FOR THE MECHANICAL BEHAVIOR OF F0AM6LAS® 


A.l Definitions 





NoDinal area over Which the load is applied, bulk dimensions. 
Load carrying cross section of cell walls. 

Nominal volume of cellular glass, bulk dimensions. 

Volume of glass in the cellular glass. 

Density of cellular glass. 

Density of glass. 

F/Ag s Bulk stress, i.e., force applied over the nominal 
dimensions. 


1 



Stress occurring in the cell walls. 

E^ Elastic modulus of the bulk material. 

Al^/Lj^ = Bulk strain produced in the material (may be 
linear or nonlinear). 


Elastic modulus of the material in the cell walls (first 
order is the value for dense glass). 




Al^/L^ s Microstrains produced in the cell wall where 
Ly characteristic dimension of the wall over which the 
strain is measured. 
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A. 2 Cellular Glass Model 


For the bulk cellular glass nstarialt 


‘b 

For the material in the cell walls: 

°w'*w ‘w 

If under load all the bulk strain is the result of extensional displaceoient in 
the wall material i.e.y no flexing or bending of the cell walls 


€ 


B 


‘W 


(A3) 


and 




(A4) 


Rearrangement leads to the following: 

’'t' " ^b/'Hi ■ f") 


Noting that lead to: 




(A6) 




Asauming an "egg crate" nodal f depicted in Figure A-1 for cellular glaaa» 
the fraction of glaaa in the noninal volune of the cellular glaaa ia the 
following I 



(A7) 


where L and i are defined in Figure Al* The fraction of load carrying glaaa 
area in the noninal area of the cellular glaaa ia given byt 





(M) 


For a unit length cell (L “ 1) Equation A7 becomea: 


i-i'-vyv. 



(A9) 


Solving for j2 yielda: 



Substituting Equation AlO into Equation A8 yielda t 


(AlO) 
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(All) 


CoBbinlng Bquaciona A6 and All yicldat 


E /B - A /A 
B W W B 





Fron Equation A6 ona findat 
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(A12) 


(A13) 


Substituting in values of O.U gm on (8.5 pcf) for the density of 
Foamglas® and 2.55 gm cm~^ (159.2 pcf) for the density of dense glass yields 
the stress of the cell wall to be 27.84 times the measured bulk stress for 
this simplified model. 
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